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1 Introduction 

The Superconducting QUantum Interference Device 
is used to detect weak biomagnetic signals since in 
the early seventies. In that time the measurement 
was based on point contact rf SQUIDs operated at 
20 M Hz, w ith to a flux noise of the order of 10" 4 
<I> J V Hz . After finding methods to fabricate reli¬ 
able thin-film devices, the dc SQUIDs started gradu¬ 
ally replace the rf SQUIDs. [1] The dc SQUID with 
an integrated signal coil had a flux noise of the order 
of 10 -5 <E> 0 /VHz. In that time the electronics was 
basedona flux modulationschemecombined withan 
cooled impedance transformer and a room tempera¬ 
ture JFET amplifier. [2] The noise of an autonomous 
SQUID was much better than that of the coupled de¬ 
vice. Owing to a thin film technology, the coupling 
constant between a dc SQUID and a signal coil has 
become high indicating that an extra noise may be re¬ 
lated to a signal coil in some way or the other. In the 
mid-eighties we introduced a method to damp the sig¬ 
nal coil resonance and to eliminate excess noise. [3] 
Mainly due to the better understanding of the SQUID 
dynamics a flux noise of tightly coupled magnetome¬ 
ters dc SQUID was reduceddown tothe level of 10 -6 
<E>o/ VHz. In the early ninentiesPTB adopted a direct 
readout based on a positive feedback method [4] and 
VTT started to use a noise cancellation method. [5] 
Neither the positive feedback nor the noise cancella¬ 
tion decrease the SQUID noise but they both make the 
electronics simplier and faster. Since a current mod¬ 
ulation is not adopted, 1/f noise from critical curent 
fluctuations cannot be eliminated without extra cir¬ 
cuitry and thus large Josephson junctions should be 
used. 

Owing to the strong progress in the eighties and in 
the early ninenties the field resolution of the SQUID 
magnetometers became good enough to meet envi¬ 
ronmental noise sources; additional field resolution 
was not necessarily needed. In the same time many 
groups started to develop high T c devices and the de¬ 
velopment of low T c dc SQUIDs was decelerated. As 
early as in the beginning of nineties we discovered 


that the damping of the Josephson junctions and the 
thus the noise can be markedly reduced. This finding 
evoked a basis for several different kinds of SQUIDs 
but they are not yet utilized in biomagnetism. [6] Re¬ 
cently, we found out that simultaneous optimization 
of the SQUID and the pick-up loop enables us to im¬ 
prove the field resolution. 

In this paper we will show that the noise of the 
SQUID based on two Josephson junction can be 
markedly improve by removing the damping of the 
junctions. In addition, we will show that the field 
resolution of the SQUID magnetometer is fundamen¬ 
tally limited by the radius of the pick-up loop ant the 
London penetration depth. The field resolution of the 
optimized device is so high that even a tiny pick-up 
loop (about 6 mm in diameter) is large enough to 
have a field resolution lower than 1 fT/Vffz. This 
all means that the future MEG instruments could have 
more channels, be more sensitive to deep sources and 
have lower price. Unfortunately, the recent demand 
is too small and the competition too weak to impel 
improvement of the multichannel systems. 

2 Noise of the SQUID based on two 
Josephson junctions 

There are several ways to provide damping for the 
SQUID circuit composed of two hysteretic Joseph¬ 
son junctions. Damping increases friction for both 
degrees of freedom in the SQUID dynamics stabi¬ 
lizing it. In Fig. (la) an arrangement where two 
impedances Z v (co) and Z (p (co), one for each degree 
of freedom, are used to damp the system. To isolate 
the current or the voltage source, impedance Z b ((o) 
is introduced. Figure (lb) shows a damping scheme 
where impedance Z^co) is set across the junctions. 
There is also an additive impedance circuit, Z dc (( 0 ), 
which is not necessarily required but enables to mod¬ 
ify IV characteristics. In general, searching for a low- 
noise SQUID is to search for minimal damping which 
ensures a stable operation. 

Following the general scheme shown in Fig. lb 
we have analyzed the SQUID design shown in Fig. 
2a. The complete analysis will be found in Ref. [7], 




Figure 1: Different possibilities to damp the double 
junction SQUID. 


The junctions are damped at high frequencies 
with an RC-circuit (R ac and C ac connected in se¬ 
ries). An ideal voltage source is connected to the 
SQUID throu gh inductance L b and external damp¬ 
ing at low frequencies is arranged by an LC-circuit 
(L dc and R c / c set in series). Assuming that (3c = 
2nR ac C 2 I c /<& 0 = 1 we can write the energy reso¬ 
lution in the form e(r) = S^/(2L) = ak B T^LC , 
where 
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Here u b = U b /(R a J c ) denotes the normalized 
bias voltage and the normalized resistance as r = 
R dc =2R ac . If r = 1, the zv-curve of the SQUID 
resembles that of the conventional dc SQUID and 
the noise a = 10. Lower values of r leads to iv 
curves with negative output resistance - the character¬ 
istics differ markedly from that of the conventional dc 
SQUID. We are interested in devices having moder¬ 
ate damping and thus we analyze here devices where 
0<r<lorr = 0. 

First we assume that r ^ 0. The normalized flux- 
to-current transfer efficiency g i(p ° = di b /dip a = 

-sin(p a cos(p ( /^u b 2 +cos 2 (p a and dynamic output 



Figure 2: Schematic diagram of the un SQUID. 
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If r ranges from 0 to 1, there is, at least, one point 
of operation where the dynamic output conductance 
becomes zero. From Eq. (2) we find that g d = 0, if 
u h = (7 - r) IcoscpJ / Vr(2-r) . Consequently, 
the voltage noise of the preamplifier has a negligible 
contribution to the SQUID noise. We call this device 
a high gain SQUID (hg SQUID). If we optimize Eq. 
(1) by taking into account Eq. (2), we find that 
a hg °P‘ - 7.6, if (p a opt = 7 1/4, r upt ~ 0.25, and $ L opt = 2. 
If we set r = 0 in Eq. (1), we find that a un opt = 
4.56 + u b 2 , if q a opt = n [l+(l-u b 2 )]/4 and b L - 
(u b /2 k) 2 . Since the SQUID has no any dc resistance, 
we called this device an unshunted SQUID. In prac¬ 
tice u b < 1, and thus ct un opt ~ 4.6. The best results 
obtained from numerical simulations lead to the same 
value for a. 

We can modify dc SQUID even more by setting 
RC-circuit (R ac /2, and 2C ac ) only across the SQUID 
output. Under this condition two last terms in Eq. 
(1) can be neglected and the optimization leads to 
q> a opt = n/2. In practise, the noise rounding pre¬ 
vents us from using this point and thus q> a opt < n/2. 
The optimization of the modified un SQUID leads to 
a un° P mod = 0.5, if $f pt = 371. In practise, (3 7 = 371 
is too high value to keep the system stable. Instead 
of using (3/ = 371 we set (3 7 = n/2, which is low 
enough to guarantee a stable operation. This choice 
leads to a u f pt mod ~\, which is near five times better 
than that of the un SQUID and more than ten times 
better than that of the conventional dc SQUID. If we 
combine the modified un SQUID and the hg SQUID, 
we get a modified hg SQUID whose noise is of the 
same order of magnitude than that of the un SQUID 



but its gain is much higher. The modified un and hg 
SQUIDs may need some high frequency damping to 
flux oscillations at the plasma resonance. This can be 
done, however, in such a way that the noise does not 
markedly increase. 

3 Noise of the Magnetometer 

In the previous chapter the noise of the autonomous 
SQUID was analyzed. In practise, the signal coil 
forms an LC-resonance and a microwave transmis¬ 
sion line leading to microwave reflections. The in¬ 
put circuit without damping increases the noise due 
to thermally activated switching between the resonant 
state induced by the signal coil and the running or 
beating mode of the SQUID. We assume that if the 
transmission line is appropriately terminated the exis¬ 
tence of the signal coil does not increase excess noise 
at all. 

Figure 3 shows the schematic diagram of the inte¬ 
grated magnetometer, where the signal coil consists 
of n turns and the pick-up loop is composed of m 
slices connected in parallel. There are many ways to 
accomplish the SQUID layout but we use this geom¬ 
etry to simplify our calculations. In addition, our new 
SQUID design based on the un SQUID (unshunted 
SQUID) resembles that shown in Fig. 3.[8] In the 
construction appropriate damping is achieved sim¬ 
ply by terminating the microstriplines formed by the 
overlaps between adjacent slices. In addition, to pre¬ 
vent reflections at the edge of the washer the char¬ 
acteristic impedance of the signal coil is m times 
smaller than that of the overlaps. 

Termination eliminates the excess noise related to mi¬ 
crowave reflections and LC resonances but it creates 
high frequency noise currents entering the SQUID 
loop. Owing to the nonlinear characteristic, Joseph- 
son junctions mix high frequency noise down to low 
frequencies. 

The magnetic field B threading the pick-up loop 
changes the flux in the SQUID ring as 

<^ = _ m 2 n 2 L^ B 2 A 2 _, (3) 

“ L 2 (1 + m 2 n 2 L g /L, + nr 2 L^, /L,) 2 

where L^is the inductance of the SQUID washer 
and L l that of the single turn pick-up loop, respec¬ 
tively. The number of turns in the signal coil is n 
and m denotes the number of sub-loops in the pick¬ 
up loop. The total stray inductance related to the mi¬ 
crostriplines formed by the signal coil andthe SQUID 



Figure 3: Layot of the SQUID magnetometer. The 
large pick-up loop (10 mm- 20 mm in diameter) con¬ 
sists of several parallel sub-loops and the signal coil 
forms the spiral microstripline on the SQUID washer 
(typically 100 pm - 200 pm in diameter). 


washer and themicrostiplines formed by the overlaps 
of the adjacent sectors is denoted Ip a r. The effective 
area sensing the magnetic field is A. 

The signal coil and the overlaps between adjacent 
sub-loops in the pick-up loop are modeled with trans¬ 
mission lines. The microwave terminations cause 
high frequency flux noise in the SQUID ring which 
is mixed down to low frequencies by the Josephson 
junctions. To eliminate reflection the termination at 
the end of each sub-loop should be Z tp = Z p = mZ s , 
where Z s is the characteristic wave impedance of the 
signal coil. 

If the SQUID ring is disturbed by high frequency flux 
<A>d then the junctions convert it to low frequencies as 
(® 2 lf )= (7t/2) 2 p 2 cos 2 (7td> 6 /O 0 ) <d> 2 , :/ ), if (3 « 
1; where (3 = 2Lf ( J<t> 0 is the normalized SQUID in¬ 
ductance. If <J> 6 = 7td> 0 /4, then <<h 2 ^) = (3 2 (<h 2 ; y). 
The mixed down noise can be given in the form 
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where L par = 2co s kln 2 l s + r^/m. Here l s and l p 
denote the inductance per unit length in the signal coil 
and the pick-up loop, r p is the radius of the pick-up 
loop, and co s is the width of the signal coil. Since 
Z s = ZJm then l p = ml s . If the microsrtipline is 
well terminated Z s = 4ljc~ = c E l s / V©, where 
0 = 1 + X L / (htf and c e the speed of the signal 
in the transmission line, \ L is the London penetration 
depth, and h is the thickness of the insulation barrier. 



The pitch of the spiral signal coil is 2(O s C2. 

It is well know that the intrinsic noise of the 
SQUID can be given in the form (^l sq ) = 
2 a yk B TL^LC ,where constant a depends on the 
SQUID type and y = (fU 2 + 4a/]3)/5. The capaci¬ 
tance of the Josephson junction is C. In the conven¬ 
tional dc SQUID optimized in respect of (3 C and (3 
parameters, a = 12. The numerical simulations and 
our analytic estimations predict a = 4-5 for the un 
SQUID. [8] Setting <D 2 = (®\ sq ) + (<i>\ pick ), we 
get 


If A ,l = 80 nm, x = 4, 4/= 100 nm, C = 1 pF , 
c 8 = 10 8 m/s, £1 = 1, and a = 12, then V<5 2 > = 
0.05(10 mm /r p ) 514 fT. Thus the magnetic field noise 
level of 1 fT, could be reached by using a pick-up coil 
having a diameter of 6 mm only. Since the distance 
from the magnetometer to sourcesin brain is typically 
about 20 mm, the optimal baseline should be much 
longer than the diameter of the loop. This means that 
an integrated planar gradiometer does not effectively 
utilize the whole chip area favoring a magnetometer. 
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where L is the SQUID inductance screened by the 
pick-up loop and ¥ = 4ak B TL 2 I A 2 . The field reso¬ 
lution related to the intrinsic SQUID noise depends 
strongly on the SQUID inductance. L can be ex¬ 
pressed as a function of the washer inductance L a as 
follows: 

L p + l s(0 m 2 n 2 

a L p + m 2 n 2 (La, + 2(o s £ll s ) ' 1 ’ 

where L p = m 2 r p l s + L t and l s(0 = 2(O s £ll s . In¬ 
serting Eq. (6) into Eq. (5) and optimizing the field 
resolution with respect to m, n and (3 we find 


- We have shown that the noise of the SQUID based on 
two Josephson junction can be much smaller than that 
of the conventional dc SQUID. In addition, the field 
resolution can be markedly improved by optimization 
the entire device as a whole; all the SQUID, signal 
coil, and pick-up loop parameters will be determined 
by the pick-loop radius. By replacing a dc SQUID 
with an un SQUID and optimizing the whole mag¬ 
netometer, a low noise device could be fabricated on 
extremely small silicon chip. Another possibility is to 
integrate several small magnetometers (e.g., five) on a 
single chip and compensate for the homogenous field 
would be by setting a separate reference magnetome¬ 
ter above the chip. An easy way to compensate for a 
homogeneous field would be to inject a feedback cur¬ 
rent of the reference SQUID through all the SQUIDs 
sensing biomagnetic signals. 


where the inductance of the circular loop is given 
as L i = 2x!i 0 Kr p , and K s is the fringe field 
factor. Equation (7) is found for $ opt = 

0.64 (V D\i 0 h C ) and h opt = = 4/. The op¬ 
timization also fixes the values for m and n. Since the 
parasitic capacitance arising from the microstipline 
does not increase the noise, the thickness of the in¬ 
sulation layer h should be as small as possible. How¬ 
ever, a thin insulation layer does not cover steps well 
enough; in practice, it should exceed the thickness of 
the lower superconducting fi lm . Thus we get a simple 
guidline h > t e j ~ X L , and the fileld resolution de¬ 
pends only on the London penetration depth and the 
radius of the pick-up loop. The optimization shows 
that all the SQUID parameters besides the junction 
capacitance are related to the radius of the pick-up 
loop and thus the magnetometer should be optimized 
as a whole. 
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